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ing volumes of cosmological size, larger than the causal horizon at the time, must 
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I. INTRODUCTION 

All current cosmological models work within the paradigm in which the present distri- 
bution of matter in the universe is the result of mostly gravitational evolution operating 
since the time of matter-radiation equality on some seed of initial and very small quantum 



2 



fluctuations in the mass (energy) density fleld of an, otherwise, homogeneous, isotropic and 
flat Friedmann-Robertson- Walker universe. On the large comoving scales of cosmological 
interest this seed of initial density perturbations at the time of equality is well described by 
the Harrison-Zeldovich spectrum of gaussian fluctuations 

The most characteristic feature of the spectrum of primordial anisotropies is the lin- 
ear dependence, (AMy)^ oc S, of the variance of mass (energy) fluctuations in restricted 
spatial sub-volumes of comoving cosmological size L on the area 5* ~ of the surface 
that bounds the sub-volume, rather than on its volume V ~ L^. This speciflc feature 
(AMy)^ oc is commonly known in the literature as scale invariance of primordial mass 
(energy) anisotropies because they produce a gravitational potential $ ~ G^^^ ~ G that 
does not depend on the comoving scale L. The scale invariance of the seed of primordial 
cosmological anisotropies at the time of equality was flrst observationally tested ten years 
ago by the COBE measurements Q of the temperature anisotropies in the cosmic microwave 
background radiation (CMBR) and has been more recently conflrmed by the high-precision 
data of the WMAP collaboration and others jj] . 

In the context of standard cosmology, nonetheless, comoving scales L of cosmological 
interest, H~^{teq) ^ L < H^^ito), are much longer, by orders of magnitude, than the 
comoving causal scale H~^{teq) at the time of matter radiation equality: H-'^{t'^\ ~ ^^"^^ 
where H~^{tQ) is the comoving size of the present causal horizon. A central problem in 
contemporary theoretical cosmology is the quantitative explanation of the mechanism that 
produced by the time of equality scale invariant anisotropies (AMy)^ ~ over cosmological 
comoving scales L, which are much larger than the causal horizon at that time. The problem 
arises from the assumption that scale invariant mass (energy) anisotropies in comoving 
volumes of cosmological size L must have been produced by fluctuations in Fourier modes of 
the density fleld with comoving wavelength of comparable size A ~ L ^ H^^{teq). We label 
this assumption as [HP] for easier reference later on. If this assumption were correct, then 
it would be justified to claim that in the context of standard cosmology we lack a causal 
explanation of the mechanism that generated, at teq, scale invariant mass anisotropies in 
comoving volumes much larger than the causal horizon at that time. This enigma is usually 
called the origin of structures problem of standard cosmology. The most outstanding of 
the mechanisms that have been invented to solve this problem is inflation, which roughly 
speaking proposes that the whole observable present universe was once causally connected in 
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the remote past and the density inhomogeneities were already imprinted at that early time 
before the universe underwent a finite period of exponential expansion which stretched it to 
its current huge size This mechanism would explain the presence at the time of equality 
of fiuctuation modes in the density field with comoving wavelength larger than the standard 
causal horizon at that instant and with the appropriate power spectrum. Obviously, the 
logic of the inflationary explanation of primordial structures also relies on assumption [HP] . 
The main purpose of this work is to bring into attention and revisit the assumption [HP] 
under which the origin of structures problem has been raised in the context of standard 
cosmology. 

Assumption [HP] is formally expressed through the fast estimation (fTTjl of integral (jlj) 
for the variance of mass (energy) anisotropics over cosmologically large comoving volumes 
(see, for example, eq. 9.1,9.2 and 9.3 in 6|). According to this estimation mass (energy) 
anisotropics over volumes V of cosmologically large comoving size L are expected to be 
generated by the dominant contribution to integral Q from random fluctuations of the 
stochastic density field (pQ) with comparably large comoving wavelength A; ~ while the 
contributions to this integral from the ultraviolet modes k ^ 1/L seemingly cancel out and 
are expected to be sub-dominant or negligible. Therefore, under assumption [HP] only a 
power spectrum with spectral index n = 1 over the range of cosmologically short comoving 
momenta k -C H{tf.g) would be able to generate scale invariant mass (energy) anisotropics 
(AMy)^ ~ over cosmologically large comoving volumes L ^ H~^(tf.q), see fll2j) . This 
statement is reported in textbooks and research papers in cosmology |2|, even though the 
fast estimation (fTT|l is since long ago known in statistical mechanics to be incorrect. The 
correct estimation of integral (j3]) over cosmologically large comoving volumes V for a power- 
law power spectrum with spectral index n is given in (fT^ . (fT^ . ()15p . for all possible values of 
the spectral index. The new estimation shows that any power-law spectrum with spectral 
index n > 1 would, in fact, generate a pattern of scale invariant cosmological mass (energy) 
anisotropics. 

Scale invariant anisotropics correspond to stochastic density fields whose random 
fiuctuations have been globally constrained to conserve the total mass (energy) of the system. 
The reason why fast estimation fails to predict the correct behaviour (fT3j) is precisely 
the fact that in these globally constrained systems the largest contribution to integral Q for 
mass (energy) anisotropics over macroscopically large volumes of comoving size L does not 
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come from modes with comparably large comoving wavelength /c ~ 1/-^, whose contribution 
is suppressed by the global constrain, but from the ultraviolet modes k > H{teq) with 
comoving wavelength within the horizon at the time of equality. The anisotropies are 
produced by a mechanism of local rearrangement of matter in the uniform universe that can 
be visualized through the following metaphor: people moving short distances (i.e. density 
fluctuations with short comoving wavelength) from town to town at different sides of a 
border can generate fluctuations in the total number of inhabitants in each of the countries 
separated by that border (i.e. mass anisotropies over cosmologically large volumes), even 
though the typical size of those countries can be much larger than the distance between the 
two towns. Obviously, the longer the common border between the two neighbour countries, 
the larger can be the size of the surface fluctuations in the number of inhabitants in each of 
the countries. Of course, the displacement of population does not change the total number of 
inhabitants in the two countries together. This mechanism is well understood in statistical 
mechanics since long ago, but for some reason it has not been clearly noticed before in 
cosmology. Statistical systems that show scale invariant anisotropies are known in the 
literature as glass-like or lattice-like. In quantum Hall effect in condensed matter physics, 
for example, low-energy excitations are known to be associated to the borders, as a direct 
consequence of the global constrain imposed by the incompressibility of the gas of electrons 
confined in two spatial dimensions Qj. 

At the light of these comments we realize that all stochastic density fields characterized 
by the same average density po and constrained to globally conserve the total mass (energy) 
of the system produce the same pattern of scale invariant mass (energy) anisotropies 
over cosmologically large volumes, even though they can be defined by very different power 
spectra. Hence, all power spectra with indexes n > 1 are, in principle, macroscopically 
indistinguishable. This degeneracy has been ignored for many years in the literature, but it 
is absolutely crucial to understand the statistical properties of primordial cosmological mass 
(energy) anisotropies at the time of matter radiation equality. 

In next sections we will revisit in detail the arguments and concepts that we have briefiy 
introduced in this section and will explicitly show that random fiuctuations in the density 
field at the same instant of matter radiation equality with comoving wavelength shorter than 
the horizon at that time H~^{tf.g) can account, when they are constrained to conserve the 
total mass (energy) in the whole system, for the scale invariant primordial mass (energy) 
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anisotropies over cosmologically large comoving volumes. 

In the context of this discussion we feel it is necessary for historical reasons to comment 
and clarify a mechanism of local rearrangement of matter at the same instant of equality first 
considered, an discarded, long ago by Y.B. Zeldovich as responsible for the mass (energy) 
anisotropies over cosmologically large comoving volumes. The apparent incapability of this 
mechanism to generate the anisotropies is cited in textbooks in cosmology proof 
of the impossibility in the standard cosmology to generate prior to the time of equality 
scale invariant anisotropies over cosmologically large comoving volumes through causally 
connected physics. The argument is also presented as a motivation for the inflationary 
explanation of the origin of primordial structures. We will show, on the contrary, that 
Zeldovich's mechanism does indeed generate scale invariant cosmological anisotropies and 
was discarded only on the ground of a misled analysis based on the wrong estimation (fTT|) . 

To model and reanalyze Zeldovich's proposal we consider a free scalar hamiltonian density, 
which is quadratic in the fundamental scalar field and its conjugate momentum. Zeldovich 
nicely noticed that random fluctuations of the fundamental fields with comoving wavelength 
within the horizon can, through the quadratic coupling, generate density fluctuations with 
comoving wavelength much larger than the horizon k^^ ^ H^^{t(.q). This is a trivial prop- 
erty of Fourier analysis. This observation, although correct, was obviously motivated by the 
mistaken assumption [HP]: the quadratic coupling offered a causal mechanism to generate 
at the same instant of equality a density power spectrum over the range of cosmologically 
short momenta and, thus, according to fast estimation (fTT|l . mass (energy) anisotropies over 
cosmologically large volumes. On the other hand, a brief and old argument that we repro- 
duce in Section V, provided also by Zeldovich, assures that the density power spectrum over 
the range of cosmologically short momenta k <^ H{t(.q) generated through the quadratic 
coupling of causally connected modes must have an spectral index n > 4 and, therefore, it 
was thought that the variance of mass (energy) anisotropies over cosmologically large vol- 
umes should decrease, at least, as (AMy)^ ~ 1/L with the comoving size L of the considered 
volume ()12j) . Hence, it was concluded that local rearrangement of matter, although able to 
generate anisotropies over cosmologically large volumes, is unable to generate scale invariant 
cosmological anisotropies (AMy)^ ~ L^. This conclusion, based on the wrong assumption 
(fTT|) . is fatally incorrect as we will show: the variance of mass (energy) anisotropies gen- 
erated through local rearrangement of matter over cosmologically large volumes are indeed 
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scale invariant (AM^)^ ~ L? (fT3j). as we should have expected because the fluctuations 
are globally constrained to conserve the total mass (energy) of the system. Zeldovich's 
analysis of his own proposal got wrong because, biased by incorrect fast estimation (fTT|) . 
it focused on the ability to produce a non-zero density power spectrum over the range of 
cosmologically short comoving momenta through the quadratic coupling in the hamiltonian 
density and missed the fact that the largest, scale invariant, contribution (AMy)^ ~ to 
cosmological anisotropics comes from random fluctuations of the density modes with the 
shortest wavelength (jlSp . while the contribution of density modes with cosmologically large 
comoving wavelength is sub-dominant, (AMy)^ ~ In other words, in this model the 

largest cosntribution to mass (energy) cosmological anisotropics does come from causally 
connected random fluctuations of the fundamental flelds when they couple to produce den- 
sity fluctuations with comoving wavelength still within the horizon, and not by modes of 
the fundamental fields that couple to produce density fluctuations with cosmologically large 
comoving wavelength. 

The irrelevance of the quadratic coupling is evident when, in the same framework set up 
to model Zeldovich's proposal, we consider the density field (0) to be proportional to the 
conjugate momentum field, instead of being given by the quadratic hamiltonian density, and 
study spatial mass (energy) over cosmologically large comoving volumes. The momentum 
field can only fiuctuate, in this example, in modes with comoving wavelength shorter than 
the horizon and, in absence of quadratic coupling, there is no way to generate density fiuc- 
tuations with wavelength longer than the horizon. Yet, an analytical evaluation of integral 
Q over cosmological volumes show that the anisotropics over volumes of comoving size 
larger than the horizon are still scale invariant. The reason ()24|) is the same as above: the 
conjugate momentum field operator is associated to a global charge of the system, which 
forces total mass (energy) conservation. In conclusion, it is not necessary to generate den- 
sity fiuctuations with cosmologically large comoving wavelength in order to produce scale 
invariant anisotropics over cosmologically large volumes. Any mechanism of local rearrange- 
ment of matter can succesfuUy generate scale invariant anisotropics over cosmologically large 
volumes. Zeldovich's old example, properly analyzed, is only a particular example. 

The analysis that we present here obviously lifts the obstacle that has long stood in the 
way to provide a causal explanation for the origin of the large scale structures of the universe 
in the context of standard cosmology (we mean standard cosmological model, without the 
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appendix of a preceding inflationary expansion) and provides an attractive alternative to 
inflation to solve the origin of structures problem. In principle, this alternative mechanism 
should probably involve only physics at the scale of matter radiation equality, Tgg ~ 1 eV, 
instead of physics at the very high energy scales up to the Planck scale that are usually 
summoned in inflationary mechanisms. Of course, this does not excluse that new physics 
could be relevant at the instant of equality to explain the origin of primordial structures. 
Moreover, the generic problem of tuning the initial cosmological conditions in inflationary 
models in order to produce a continued period of exponential expansion simply vanishes 
in the context of the alternative scenario that we suggest: in this alternative scenario it is 
demanded in a naturally simple way that the fluctuating density field is in its ground state 
at the instant of equality. 

The paper is organized in nine sections. In Section II we review the concepts and tools 
needed to describe the statistical properties of density anisotropics at the time of matter 
radiation equality in the homogeneous, isotropic and flat FRW universe. In Section III we 
characterize the power spectra of gaussian fluctuations of the density field that produce scale 
invariant mass (energy) anisotropics over cosmologically large comoving volumes. It is ex- 
plicitly shown that any power-law spectrum V{k) ~ Ak"' with spectral index n > 1 produces 
an scale invariant spectrum of fluctuations (AMy)^ ~ in volumes of cosmologically large 
comoving size, L > H~^(t(.fj). This claim is the main reason why the analysis performed in 
this paper disagrees with the main conclusions that are found in the literature, where it is 
claimed that only a power-law spectrum with spectral index close to n = 1 can produce scale 
invariant mass (energy) anisotropics in cosmologically large comoving volumes. The motives 
of our disagreement are extensively discussed in this section III. The calculations presented 
also prove the central claim of this paper, that random fluctuations at the same instant of 
equality in the Fourier modes of the density field with comoving wavelength shorter than 
the horizon at that time can account for the observed scale invariant primordial anisotropics 
over cosmologically large comoving volumes. We go on in Section IV to discuss the charac- 
terization of scale invariant mass (energy) random anisotropics in terms of the two points 
correlation function of the fluctuating density field. We can advance here that we find that 
scale invariant anisotropics correspond to a certain class of short range correlation functions 
that result when the global mass (energy) contained in the system is conserved, and there- 
fore, not allowed to fluctuate. The aim of section V is mainly to clarify the arguments that 
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historically laid the so-called origin of structures problem [HP] of the standard cosmology. 
We discuss these arguments in the setup of a quantum field theory and describe the mech- 
anism of generation of scale invariant primordial cosmological mass (energy) anisotropies 
by local rearrangement of matter as random quantum fluctuations of the density field in 
the homogeneous and isotropic ground state of the system at the same instant of equal- 
ity. Notwithstanding these first five sections of the paper constitute a self-contained body, 
which justifies and proves the claims that appear in the abstract, introduction (Section I) 
and conclusions (Section IX) of this paper, we have judged appropriated to include three 
additional sections where we introduce concepts that shall be relevant to set up a complete 
quantum theory of density anisotropies coupled to metric perturbations of the homogeneous 
and isotropic FRW space-time background. In Section VI we discuss technical concepts on 
the renormalization of ultraviolet divergent expressions in the theoretical setup to conclude 
that, in general, the size of the scale invariant spatial anisotropies can only be defined as 
an external renormalizable parameter of the theory, like masses or coupling constants. In 
Section VIII we obtain the renormalization equations that describe the dependence of this 
parameter on the resolution scale at which the system is probed. In Section VII we review, at 
the light of our findings, the linearized gauge invariant formalism |6] of density anisotropies 
coupled to metric perturbations of the FRW expanding space-time background. Section IX 
summarizes the conclusions of this paper. 



II. DENSITY ANISOTROPIES IN HOMOGENEOUS AND ISOTROPIC 
STATISTICAL SYSTEMS: CONCEPTS AND DEFINITIONS. 

We review in this section some of the concepts and tools of statistical mechanics that 
are widely used to describe the cosmological density anisotropies in the early universe at 
the time teq of matter radiation equality. Our aim is to gather some important results from 
statistical mechanics upon which we will develop our discussion. This review will also help 
us to fix the notation that we use in the rest of the paper. For the sake of simplicity we 
work in comoving coordinates and fix the scale factor due to the expansion of the universe 
to be equal to one at the time of equality, a(teq) = 1. 

The statistical system is described by an homogeneous and isotropic density field p{x) in 
3D flat space n = R^: 
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p(^) = po + poJ^, k (1) 

If the fluctuations are small the coefficients 5^ = 5*^ are independent random complex 
variables with normal distribution, {5*^ 5j:^) = {27c)^V{ki)6^{ki — ^2), and zero expectation 
value, (5^) = 0, so that each local variable p{x) is real and has expectation value (p(x)) = po 
independent of the position x. This average value is assumed to be non zero, po 7^ 0. The 
central limit theorem implies that the fluctuations of the density field p(x) are gaussian. 
The function V{k) > 0, which gives the variance of the fluctuations of each of the random 
Fourier modes is usually called the power spectrum of the statistical fluctuations and it 
is the most common statistical tool used to describe cosmological models. If the system is 
isotropic, the power spectrum V{k) = V{k) depends only on the modulus of the momentum 
that labels each mode, k = \k\. 

Statistical fluctuations of the density field in macroscopic, but finite, sub- volumes V C 
are usually described in terms of the statistical magnitude cr'^iV) = -^((^/y rf^x ^ ). 

This magnitude is sometimes called the squared density contrast over the volume V and 
de„oted ,y {^%. I. ...s of ... powe. spec... of *e fluc.uaf.o,. .he s,ua..ed de,..y 
contrast can be expressed as: 

where the geometric factor Fv{k) is given by the expression 

Fv{k) = f d^xe-'^-\ (3) 

J 

We wish to notice that this geometric factor Fy{k) is the restricted integral over the sub- 
volume V of the Fourier modes of expansion (0) in flat 3D space. 
Step by step the proof of relation proceeds as follows: 



V^\Jv Po / Jv Jv J (27r)3 J (27r)3 ^ / 



The statistical variable M{V) = Jy d^x p{x) describes the total mass (energy) con- 
tained in the sub-volume V. This is the mathematical object that describes the macro- 
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scopic properties of the statistical anisotropies in the density field (H)). Its average value 
{M{V)) = PqV is proportional to the volume V of the considered region. Its variance 
(AMv)' ^ {[MiV)-{MiVm = {M\V)) - {M{V)f = {[Jy d^x (pix) - po)]') = 
cr^(V^) X (M(V^))^ measures the typical size of its gaussian fiuctuations. Combining the 
previous expressions it is straightforward to obtain: 

{AMyr = plJ-0^^V{k)\Fyik)\\ (4) 

The last tool we want to introduce in this section is the two points correlation function 
of the density field, F{x, y) = (p(x) p{y)) — pi- If the system is homogeneous the two points 
function must depend only on their relative position r = x — y. If the system is also isotropic 
then the two points function F{r) can depend only on the distance between the two points. 
In terms of the power spectrum, we can write 

m =pII^s V{k)e^^^-\ (5) 

The proof is straightforward: 



pI I ^ / d% V{k,)6'\h - k2)e+'''-'e-^''-' = pU Ms V{k)e+^'<^-y^ 



(27r)3 J - • " (2^)3 

The variance of mass fiuctuations (j3)) can then be conveniently expressed in terms of the 
two point function: 

(AMy)2 = 1^ d^x d'y Fix,y}. (6) 

Let us remark at this point that V{k) must be integrable over the whole 3D momentum 
space. 



d^k r{k) < oo, (7) 

if the two points function ^ is well defined at r = and the squared density contrast 
vanishes when integrated over very large volumes, limv-, q = 0. This condition (0) 
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allows the spectrum to diverge at = 0, with the only condition that, if it diverges, it does 
slower than 

= 0. (8) 

The discussion in this paper will focus in the statistical properties, namely power spectra 
and two points correlation functions, of scale invariant systems, in which the variance of mass 
(energy) fluctuations in restricted sub-volumes V is proportional to the area of the surface 
S that bounds the sub-volume, (AMy)^ oc S, rather than to its volume V. If the considered 
sub-volume is, in particular, an sphere of radius L those statistical systems produce mass 
fluctuations 



limk^o V{k) 



{AMvf oc Sr^L^. (9) 

Fluctuations in the density field p{x) seed in turn fluctuations in the gravitational po- 
tential through the Poisson's equation V^$(x) = 47rG'(p(x) — po)- When the mass (energy) 
fluctuations in cosmologically large comoving volumes are scale invariant (jH)) the fluctuations 
in the gravitational potential 

A^v- oc ~ G, (10) 

do not depend on the scale L. Statistical systems with this remarkable feature are named 
after Harrison and Zeldovich, who first identified them 



III. POWER SPECTRUM OF SCALE INVARIANT DENSITY ANISOTROPIES. 

Our aim in this section is to characterize in terms of their power spectrum V{k) the 
statistical systems whose mass (energy) fluctuations over cosmologically large comoving 
volumes have scale invariant variance ©• The two magnitudes are directly related by 
equation (jl)). 

A fast, but incorrect, estimation of relation Q which, nevertheless, appears in many 
references in the literature (see, for example, eq. 9.1,9.2 and 9.3 in j^) gives: 

r d^k ^ 1 

{AMyf =PIJ7^, V{k)\Fy{k)\' ~ pI- V{k ~ l/L)|Fy(0)p = 
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= pI^, V{k ~ 1/L)V' = pI V{k ~ 1/L) V, (11) 

where L is the hnear size of the considered volume of integration. This estimation is derived 
after noticing that the geometric factor Fv{k) evaluated at = measures the volume 

~ of the considered spatial region, Fv{k = 0) = V, and it decreases fast to zero 
for values of the momentum k larger than the inverse 1/L of the typical linear size of the 
considered volume, Fv{k > 1/L) ~ 0. 

According to this estimation a power-law spectrum with spectral index n, V{k) ~ Ak^, 
over the range of comoving momenta k < H{teq) produces mass (energy) fluctuations 

{AMyf ~ pI V{k ~ 1/L) L'r^pl^L' = plA L^"", (12) 

in cosmologically large comoving volumes, of size L > H~^{tf.q). Therefore, a power-law 
spectrum with spectral index n = 1, V{k) ~ A fc, over the whole range of comoving momenta 
k < H(teq) seems to be a necessary and sufficient condition to generate a pattern (jH)) of scale 
invariant mass (energy) anisotropics over cosmologically large comoving volumes. 

Notwithstanding this result (fT^ is cited in textbooks and reviews on the subject 0, llj, 
we choose to check it, both analytically and numerically, integrating equation (Q. It will 
turn out that p2|) gives a reliable estimation only when the spectral index n is in the range, 
—3 < n < 1, but the estimation gets fatally wrong for larger values, n > 1. Notice, in 
particular, that (fT^ predicts an scaling law (AMy)^ ~ L", with a < 2 for values of the 
spectral index larger than one, n > 1. For example, for n = 4 it predicts (AMy)^ ~ 
Clearly, this prediction must be incorrect because it has been rigorously proved that a > 2 
the scaling power a is necessarily equal or larger than 2 for any homogeneous and 
isotropic statistical system. The fast estimation (fT^ misses that for values of the spectral 
index n > 1, in spite of the rapid decrease to zero of the geometric factor \Fy{k)\'^ for values 
of comoving momenta k ^ 1/L, integral is still dominated by the contribution of the very 
large momenta k ^ 1/L, rather than by the modes k ~ 1/L. This ultraviolet divergence 
has been noticed before in the literature on cosmological density anisotropics |l^. Once this 
contribution is taken into account it is found the correct behaviour of the variance of mass 
fluctuations, 



If nG(-3,0], (AMy)2oc yi-"/^ _ (^3^ 
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If ne(0,+l], (AMy)2 OC 53/2-n/2 _ (^^^ 

If ne[+l,oo), (AMy)2oc S ^ L"^ . (15) 

These correct results ()13I14I15|) . obviously, respect the bound a > 2. We wish to remark 
that these results are not new neither unexpected. They were previously reported, for 
example, in [3] in the context of a discussion of the spectrum of primordial cosmological 
anisotropics. They are also very well-known in condensed matter physics in the description 
of glass-like systems. We rederive them here only with the aim of convincing the reader 
about their validity. As we will show, the disagreement between the old estimation (fT^ and 
the correct result (fT3j) listed above is the only reason why the conclusions of this paper do 
not agree with those widely accepted and reported in the literature. 

In order to justify the result ()15|) we now consider as volume V of integration an sphere 
of comoving radius L. The rotationally symmetric geometry of the considered spherical 
sub-volume does not change the qualitative features of the results that we want to prove, 
but it allows to obtain the geometric factor (jH]) analytically: 

Fv{k) = AtcL^-^- {sin{kL) - (kL) cos{kL)) . (16) 
[kL)-^ 

When this expression is introduced in Q we obtain the relation: 

{AMvf = 8pl J dikL)^^^^ {sin{kL) - (kL) cos{kL)f . (17) 

It is important to notice that X(w) = {sin{w) — w cos{w) f, where w = kL, is not singular 
at w = as the factor could induce to think, because {sin{w) — w cos{w)Y = 0{w^) 
when w — > 0. 

For the power spectrum we now assume a power-law dependence V{k) = -A-^^ over a 
large, but finite, range of comoving momenta < k < k^, and beyond which the power 
spectrum is cutoff by some supression factor. Later on we will specify the scale kc of the 
cutoff procedure and will characterize it as a regularization step of an ultraviolet divergent 
expression in the context of the renormalization programme of a physical parameter. For 
the moment it is enough to say that this scale can be somewhat of the order of the Hubble 
horizon at the time of equality, kc ^ H{teq). Notice also that coefficient A is dimensionless 
and can be fixed independently of the spectral index n. For the sake of simplicity we will 
assume in the analytical proof that we present below that the cutoff at kc in the power 
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spectrum is sharp, > kc) = 0. More realistic exponential or polynomial cutoffs are 
evaluated numerically. The results are presented in Fig. 1. They show the same behaviour 
fllHpi4|15|l that we observe from the analytical discussion of the simpler case with a sharp 
cutoff. 

We now discuss separately the four different cases when the spectral index G i? in the 
power-law is: case I) n > 1; case II) n = 1; case III) < n < 1; and, finally, case IV) 
—3 < 77, < 0. The spectral index n G -R is restricted by condition (jHJ to be larger than 
n > -3. 

In case I), n > 1, the variance (|T7jl can be estimated with very good accuracy by the 
analytical expression 

i^Myf-^PlL'£^'d{kL)^^, (18) 

if L ^ k~^. The reader can visualize the condition for this estimation after noticing that 
the factor ^^^^2 {sin{kL) — {kL) cos{kL)Y in the integrand of ()17|1 approximately halves the 
area left under the curve When the last integral is performed we obtain 

M ^ ,d L» ^ ,4 ^^1,L^ . ,1 A'l^L^ (19) 

where A' is a dimensionless factor that fixes the absolute amplitude of the fluctuations. 
Notice the promised linear dependence of the variance on the area S* oc of the surface 
that bounds the sphere of integration. It is also interesting to realize that integral ()17|) is 
convergent at kL —>■ 0, but it is divergent in the modes with very large momenta k ^ k^^ 
1/L close to the cutoff. The explicit way how the cutoff in V{k) regularize this expression 
is not qualitatively relevant in this discussion, as it is shown in Fig. 1, but the need for 
regularization proves that the largest contribution to the variance comes from fluctuations 
of the Fourier modes 6k of the density field (P) in the ultraviolet range k ~ kc- We also 
can understand this result by noticing that, according to (fT^ . the contribution to integral 
(0} over comoving volumes of cosmological size L coming from modes with comparable 
wavelength k 1/L decreases for n > 1 as ~ L'^"", which is subdominant to the contribution 
fll5|) ~ L^ from the ultraviolet modes close to the cutoff k ~ fec- 
it is also important to notice how the ultraviolet cutoff k^ enters expression ^9^ : this 
regularization cutoff only modifies the absolute size of the mass (energy) fluctuations, but 
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it does not modify the scale invariant linear dependence of their variance on the area 5* of 
the surface that bounds the considered sub- volume (P). We will discuss this point in more 
detail in next Section VI and Section VIII. 

Before going on to cases II, III) and IV) let us further clarify the result we have 
just proved. We have shown that fluctuations in the modes 5k of the density field (Q) with 
the shortest comoving wavelength k > H{teq), can generate scale invariant mass (energy) 
anisotropies over much larger volumes of comoving cosmological size L ^ H~^{t(,q). We have 
been asked where is the 'non-linear' mechanism that leads to this result ? The answer is in 
the integral expression for the variance of mass anisotropies: the reader can clearly see 
that fluctuations in all the modes 6k over the whole range < k < kc add up their positive 
contribution dk k"^ V{k) \Fv{k)\^ > to the total variance. It is commonly assumed that 
for a given volume V of cosmological comoving size L ^ the overwhelmingly largest 
contribution to this integral Q comes from the modes k ~ L^^ (see, for example, eq. 
(9.3) in j^). This is the assumption [HP], which leads to the 'linear' relationship (|TT|l 
between the variance of mass (energy) anisotropies in a volume of comoving size L and 
fluctuations V{k ~ L^^) of the modes 5k with comparable comoving wavelength L ~ k^^. 
We have explicitly shown above that this 'linear' relationship does not hold in case I), when 
fluctuations of the the ultraviolet modes /c ~ fcc produce the largest contribution to the 
global integral Q and generate scale invariant mass (energy) anisotropies. The 'non-linear' 
connection between mass (energy) anisotropies in cosmologically large comoving volumes 
and fluctuations in the Fourier modes of the density field with much shorter comoving 
wavelength is intrinsically wired into the integral expression (Q, which stands as the core 
hardware of the formal description of the macroscopic properties of gaussian fluctuations 
in the density field (^. This 'non-linear' connection is relevant in case I), when spatial 
anisotropies are scale invariant, but has been carelessly disregarded in the literature on the 
origin of cosmological primordial density anisotropies. 

In case II), n = 1, the integral (fT^ is apparently divergent not only at A;L ^ 1 but also 
at kL —* 0. Now it is important to remember that, as V{0) = for any n > 0, the integrand 
V{k)I{kL) in (fTTj) . of which (fTHj) is only an estimate, is regular at kL —>■ 0. The variance 
is still dominated by the ultraviolet modes k k^ ^ 1/L, but now they contribute only a 
sub-leading logarithmic divergence. 
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{AMvY ~ ^pIA — ln{KL). (20) 

To leading order the variance (AMy)^ ~ shows the same hnear dependence on the area 
of the boundary surface that was obtained for the variance in case I). 

In case III), < n < 1, the integrand in ()17j) decreases fast enough in the large k modes 
to render the integral over these modes convergent. A proper estimation shows that now the 
result agrees with the estimation in (fT^ . (AMy)^ ~ p^B' -t^L^~^, because the integral 
now is dominated by the 'linear' modes kcL ~ 1. 

Finally, we consider the case of negative values of the spectral index, —3 < n < 0. In this 
case, the integrand in the l.h.s of (jl7p decreases even faster to zero in the ultraviolet regime 
k ^ kc ^ 1/ L and the integral is dominated by the modes fcL ~ 1. The fast estimation 
presented in (fT^ is, therefore, also valid in this case. In fact we obtain, 

(AMy)2 ~ pIB"^^L'-. (21) 

This completes our proof of (|13ll4ll5p . These results show that any power-law spectrum 
over a finite range of comoving momenta < k < kc, that vanishes at k = 0, i.e. V{0) = 0, 
and whose first derivative at the origin also vanishes, i.e. V'{0) = 0, produces a pattern of 
scale invariant mass (energy) fluctuations in any volume whose comoving size is much larger 
than the comoving length scale of the cutoff, k^L ^ 1: 

{AMvY <xS r^L'^ ^ (P(0) =0 V V'{0) = 0) . (22) 

If only the power spectrum vanishes at zero, 'P(O) = 0, but the first derivative does 
not, V'{0) 7^ 0, we still have surface dependance (fT^ for the variance of mass (energy) 
anisotropics, (AMy)^ oc S"', although the polynomial growth is faster than linear 1 < 
7 < 1.5. On the other hand any power spectrum which does not vanish at A; = 0, i.e 
V{0) 7^ 0, produces in spatial regions whose typical comoving size is much larger than 
the comoving length scale of the cutoff kcL ^ 1 a pattern of mass (energy) anisotropics 
(AMy)^ oc V^, with 1 < (3 < 2, proportional to some power of the volume. We will 
further comment on this point in next section. Now we do not want to derail from the main 
discussion of this section that shall focus on power-law spectra with spectral index n > 1 or 
any linear combination of them, which produce strict scale invariant anisotropics over 
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cosmologically large comoving volumes through the dominant contribution of the Fourier 
modes with comoving wavelength shorter than the horizon k > H{tf.g). 

The result that we report in equation ()15|1 has further implications that we want to em- 
phasize: power spectra with different spectral indices n > 1 produce gaussian fluctuations 
in the macroscopic variable M{V) = Jy p{x) with the same average value poV and the same 
characteristic variance (AMy)^ = p^A!k~'^L'^ in any macroscopic volume V , for their ampli- 
tude can independently be flxed through the dimensionless prefactor A! . Therefore, we must 
conclude that power spectra with different spectral indices n > 1, or linear combinations of 
them, are macroscopically indistinguishable in the scale invariant mass (energy) anisotropies 
AMy that they produce in cosmologically large spatial sub-volumes. 

If all different power spectra V{k) ~ A;", with index n > 1, produce mass (energy) 
anisotropies in macroscopic volumes with the same variance, which scales linearly with the 
boundary area (jHl), all of them produce scale invariant gravitational potentials (fTUj). Hence, 
there is no way to formally distinguish between power-law spectra with different n > 1 
through the macroscopic mass fluctuations they produce, neither through their gravitational 
effect on large macroscopic distances. 

Does this discussion mean that power spectra with different spectral index n > 1 are 
physically indistinguishable ? Not exactly, in principle. Different power spectra would 
produce different two points correlation functions -F(r), but the differences are hidden when 
the macroscopic magnitudes (AMy)^, given by ©, are compared. We will come back to 
this point later in Section VI and Section VIII. 

We flnish this section with an estimation of the density contrast {5p/p)L = \f^L 
spherical volumes of comoving size L associated to scale invariant mass (energy) anisotropies 
(Uni). We have: 

2 {AMyf pI A'^L' A' 1 



(My)2 p2^2 (4^/3)2 ^k^iy 

(47r/3) (fccL)2 ^■^'^ ^ (kaLf 



Therefore, {dp/p)^ ~ r/^i (h \ -)2 ~ 0-25 V^tt-tw- Present bounds on the density con- 



trast over cosmological scales constrain (Sp/p) ^ 10~^. Therefore, we can constrain the 
scale kc of momentum cutoff to be such that k^L > 100 (A'Y^^ for all cosmologically large 
comoving scales, H{teq)L ^ 1. Both expressions together demand that A! < {kc/H{teq))'^, 
so that kc can be naturally larger than H{teq). 
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IV. SCALE INVARIANT DENSITY ANISOTROPIES AND TOTAL MASS 

CONSERVATION. 

We want to obtain in this section the conditions that characterize the two points correla- 
tion function (0) of statistical systems which show scale invariant mass (energy) fluctuations 
Q over macroscopically large, but finite, spatial sub-volumes V G Q of the whole system. 
These conditions will be then used to prove that random fluctuations of the density field (QJ 
produce scale invariant mass (energy) anisotropies, (AMy)^ oc S, in cosmologically large 
comoving volumes if an additional constrain of total mass (energy) conservation is imposed: 

{AMnf = {AMv^ocS. (24) 

Let us start then characterizing the two points function of scale invariant density fluctu- 
ations. Equation (jSj) can be inverted and written as 

p(^) = L f d^re-'^-^'Fif) = % Tdr r^'-^!^F{r). (25) 
Po Po Jo kr 

The power spectrum is by definition a positive function, V{k) > 0. This condition 
constrains any possible choice for the two points correlation function F{r). Yet the most 
interesting constrain is derived from the condition that appears in the right hand side of the 
logical equivalence (j^ : V{k)k=o = and ^^^|fc=o = 0. 

First, we notice 

1 

V{0) = ^ dr r^F{r) = 0. (26) 

Po -^0 

This is the only non-trivial constrain that we need to request on the two points correlation 
function F{r) in order to get scale invariance mass (energy) anisotropies according to (j22j) . 
because the condition on the first derivative at the origin k = is trivially fulfilled if the 
linear operators ^ and dr commute when acting on the function g(k, r) = r^^^^^^F(r): 

P5 dk \Jo kr J Jo dk \ kr J 

(27) 

We will not discuss here the formal mathematical requirements needed to get this technical 
condition (P7j) satisfied and will assume, for simplicity, that it is trivially fulfilled in the 
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setup we have defined. In fact, in all the standard definitions introduced in Section II it 
was assumed that linear operators in real r-space and momentum fc-space do commute. Of 
course, a further consideration of these formal aspects would be interesting. 

The logical equivalence (j22I)j therefore, implies that condition J^dfF{r) = in the 
r.h.s. of completely characterizes the two points function F{r) of statistical systems 
which produce scale invariant mass (energy) anisotropics. This condition on the two points 
function is known in statistical mechanics to characterize glass-like systems. In quantum 
Hall effect in condensed matter physics a similar constrain on low-energy excitations is 
associated to the incompressibility of the gas of electrons confined in 2D spatial dimensions. 
These excitations are known to be associated to the borders j^. 

The condition above can be written in a different way as 



[ drF{r) = ^ [ dr{{p{x) p{x + r)) - {p{x)){p{x + f))) = 
Pq Jn Pq Jn 

= \{{p{x)Mm)-{p{^)){Mm) = o, 

Po 

where M{fl) = dxp{x) is the total mass in the whole space Vt. The condition of scale 
invariance (j2SI); therefore, is satisfied if and only if the total mass (energy) M[VL) is not 
correlated to the specific value of the density field at any point x: 

{p{x)Mm = {p{x)){M{n)). (28) 

Furthermore, if we integrate the previous expression over a finite volume V we obtain that 
in systems that show scale invariant random mass (energy) anisotropics 0, 

{M{V)M{n)) = (M(\/))(M(^])), (29) 

the total mass (energy) contained in the system, M{Q), is not correlated to the mass (energy) 
contained in any finite sub- volume V G Q inside it. Of course this condition is easily fulfilled 
if the total mass M{Q) is constrained to a constant value and cannot fluctuate. Hence, 
we conclude that, in the context of standard cosmology (again, without inflation) random 
fluctuations of the density field Q are necessarily scale invariant Q if the total mass of the 
system is constrained and it is not allowed to fluctuate This statement is expressed 
through the logical implication ()24|). 
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The surface dependence of the variance of mass (energy) fluctuations in a system whose 
global energy does not fluctuate is easily understandable: mass (energy) fluctuations in 
any two complementary closed sub- volumes V and V = cl{Q — V) must be correlated, 
such that {AMvf = (AM^)2, because Mq = My + My is an extensive magnitude. This 
correlation implies that the variance of mass fluctuations in sub-volume V must be only a 
function (AMy)^ = f{S) of the area of its common boundary with V. This is what we 
found in (fT^ and (fT3j) : when n > global mass (energy) condition (j^^ is satisfled because 
V{0) =0. Furthermore, technical condition ^I7\i guarantees the hnear dependence (fT3j) once 
this condition is fulflUed. 

V. COSMOLOGICAL DENSITY ANISOTROPIES FROM SHORT 
WAVELENGTH REARRANGEMENTS IN THE DISTRIBUTION OF MATTER. 

At this point of our discussion we feel necessary to revisit and comment an old argument 
flrst discussed by Y.B. Zeldovich joj and, after him, by many other authors in textbooks j?] 
and research papers [l^. The argument discusses, under assumption [HP], the possibility 
that causally connected processes at the instant of matter radiation equality could account 
in the standard cosmology (again, without inflation) for the cosmological mass (energy) 
anisotropies at that time. 

According to mistaken estimation PHI . derived from assumption [HP], the mechanism 
responsible for scale invariant mass (energy) anisotropies over cosmologically large comov- 
ing volumes at the instant of equality would need to generate a density power spectrum 
V{k) ~ k over the range of cosmologically short comoving momenta. That is, the mech- 
anism would need to justify: flrst, the presence at the instant of equality of fluctuating 
modes in the density fleld with comoving wavelength much larger than the causal horizon 
at that time; second, the linear growth with the momentum scale of the power stored in 
these modes. The flrst obstacle seems easily superable: trivial Fourier analysis is enough 
to notice that fluctuations at the same instant of equality in modes of weakly interacting 
flelds with comoving wavelength shorter than the horizon can couple through the quadratic 
kinetic terms of the hamiltonian density to produce density fluctuations with cosmologically 
large comoving wavelength. 

On the other hand, the second requirement seems unattainable. The power spectrum 
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that quadratic coupling creates over the range of cosmologically short comoving momenta, 
k <^ H{teq), is not hnear, but quartic, in the momentum of the Fourier modes: V{k) ~ k^. 
The exphcit calculation can be found in the textbooks. We reproduce it below in this section. 
For now, the following fast estimation will suffice: 



If the total mass is conserved, then the first, /c-independent, term of the expansion cancels 
out, which implies 6k ~ 0{k) that the power spectrum of the fiuctuations V{k) = l^^p ~ 
0{k'^) and, therefore, V{0) = 0. This is in full agreement with our results (|24|) and (j22|) of 
the previous section. In addition it is also noticed in the literature that the linear term in k, 
second in the above expansion, must also vanishes if the center of mass of the distribution 
does not fiuctuate. In other words, conservation laws of the mass distribution demand 
6k ~ 0{k'^). In consequence, the power spectrum of mass (energy) fiuctuations generated 
locally grows, at least, as the quartic power n = 4 of the momentum, V{k) = ~ A;^, 
over the whole range of cosmologically short momenta, k < H{teq). 

Introducing this quartic power spectrum in estimation p2|l . derived from assumption 
[HP] through eq. we obtain the clear prediction (AMy)^ ~ that the variance of 

mass (energy) cosmological anisotropies generated by local modes of a fundamental field that 
couple quadratically in the density field seem to decrease at least inversely proportional to the 
comoving size of the considered volume. On the ground of this analysis it is usually claimed 
in cosmology the impossibility to causally generate in the standard cosmology scale invariant 
mass (energy) anisotropies (AMy)^ ~ L"^ over cosmologically large comoving volumes. This 
conclusion is known as the origin of structures problem of standard cosmology. We wish 
the reader to realize how heavily this conclussion relies on estimation ()11| ) .()12 |) . which was 
derived from assumption [HP]. 

Let us now show the inconsistency of the previous analysis, highlighting the use of as- 
sumption [HP] and the assessment (j3U|) for the power spectrum in Zeldovich's proposal, 
V{k) ~ k'^ over the range of cosmologically short momenta. We have learned in Section HI, 
case I, that estimation [HP], written down in eq. (fT^ . for mass (energy) anisotropies with a 
quartic, n = 4, power spectrum V{k) ~ k*, gets absolutely wrong because it does not prop- 




(30) 
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erly take into account the contribution from density fluctuations with very short comoving 
wavelength. In fact, according to the correct estimation (fT3jl the anisotropies generated by 
the power spectrum with spectral index n = 4 are indeed scale invariant, as we should 
have expected since the very beginning of this section from the global constrain on total 
mass (energy) conservation (j2H). 

Zeldovich's proposal is only a particular example of a more general mechanism of gener- 
ating scale invariant mass (energy) anisotropies over cosmologically large comoving volumes 
through the causally connected local rearrangement of matter at the time of equality in 
a previously uniform universe. How does this mechanism work? Displacement of masses 
through the border that bounds the considered sub-volume changes the total masses con- 
tained both inside and outside that volume, although the total mass (the sum of the mass 
inside plus the mass outside) is conserved and, therefore, the cosmological mass anisotropies 
it generates must be scale invariant ()24|) . 

This mechanism works even when the density fleld is considered to be linear, instead 
of quadratic, in the fundamental flelds, as we explicitly show below in this section. In 
fact, the capability in Zeldovich's original proposal to generate density fluctuations with 
cosmologically large comoving wavelength through the quadratic coupling of fundamental 
flelds is only a secondary feature without relevance, because the contribution of these very 
long modes to integral (0)) for the variance of mass (energy) cosmological anisotropies is 
negligible compared to the contribution of modes of the density fleld with short comoving 
wavelength within the horizon. 

In order to make these arguments explicit and absolutely clear we want to show how Zel- 
dovich's mechanism works in a toy-model and point out where the analysis of the mechanism 
gets wrong in the literature, under assumption We consider a massless and free quan- 
tum scalar fleld (f){x^) in the homogeneous, isotropic and flat FRW space-time background 



where x'^ are comoving coordinates and a{xo) is the scale factor of the expansion, which we 
flxed to be equal to one at the instant of equality, a(Xgg) = 1. Let us state clearly that 
the scalar fleld we use here is not the inflaton, neither it is any new particle that we are 
inventing. Action (|^T| is simply a toy-model of the mechanism proposed and discussed by 





(31) 
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Y.B. Zeldovich long ago. We could have chosen to show this mechanism in a system of 
fermions or gauge bosons. 

We are interested in describing the spectrum of energy (mass) anisotropics over cosmo- 
logically large comoving volumes due to local displacements of matter within the causal 
horizon at equality, which we describe as quantum fluctuations in the Fourier modes of the 
scalar field 



IkjR (27r)3 

with comoving wavelength shorter than the causal horizon, < H~^{teq). Therefore, in 
our toy-model we introduce an infrared cutoff in momentum km ~ H{teq) in the phase space 
of the scalar quantum field. The way how this infrared cutoff is introduced is not relevant, 
it can be a sharp cutoff or a mass term in the action (jM]) of the scalar field. The important 
point is that we want to center the attention on the role of the under-horizon modes at the 
time of equality in the generation of mass (energy) anisotropics over cosmologically large 
comoving volumes. In Section VII this toy-model is embedded into a physically motivated 



d.cnp«on of ..e a..ot.opie. t.e .au.e .v.ia,. fo.n,a,.. fl. The. we wiU 

comment on the role of super-horizon modes of the fluctuating field ()32|) . 

The dynamics of Fourier modes of (jH^ with comoving wavelength within the horizon is 
not affected by the expansion of the universe and can be described by a free hamiltonian in 
fiat Minkowski space-time, 

U = I d^x h{x) = f d^x (po + : ti*{x)ti{x) + dicf)* {x)di(t){x) :) , (33) 
where 7r*(x) is the conjugate momentum of the scalar field: 

= / ^sMPIe-^'-' (34) 

and Po 7^ is the average energy density. Only for the sake of elegance we have separated 
the average energy density as a free parameter in (jHSI) and canceled out the expectation 

value of the second term, enclosed between signs : : of normal ordering. 

The conjugate operators and obey canonical commutation relations [Qj:^,P!j = 
i {27T)^6'^{ki — ^2). Introducing the Fourier expansions (jH^ and (jMj) in (j^Hj) we can describe 
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: of 



the dynamics of the scalar field as a collection Ti = Ho + : /^^^ ^2;;^ 1^1 -^|-^a- + ^*k^k 
free harmonic oscillators. The zero mode is TYq = {Jnd^^ Po)- 

If the expansion of the space-time background preceding the time of matter radiation 
equality is perfectly adiabatic, the scalar quantum field reaches that instant in its funda- 
mental state |0). In other words, we assume here that there is not cosmological produc- 
tion of real 0-particles. The ground state of the scalar field is the tensorial product of 
the vacua of the decoupled harmonic modes. Fourier modes do fluctuate quantum me- 
chanically in their ground state -P^|0) = i <5^|0) with normal distribution of covariance 
{Q*j: Qk^ — (27r)'^5^(A;i — and can generate gaussian fluctuations in the macroscopic re- 
stricted observables Tiy = Jy d^x 7i(x) that describe the energy (mass) contained in macro- 
scopic spatial sub-volumes V, 



because the ground state |0) of the quantum field is not an eigenstate of these macroscopic 
bulk operators. 

The typical size of the macroscopic fluctuations is estimated by the variance (AEy^"^ 



I' 



(0|7iy|0) — (Oj7iy|0)^. A formal expansion of this expression was obtained in and 
rederived in (l7l |. 

{AEvf = ——. / d'hd'k2\Fv{h-h)\'\h\\k2\{cos{e)-l)\ (35) 
4(27rj'^ JkjR 

where kuv ~ H{t(.q) is some cutoff that regularizes the theory in the ultraviolet and 

cos{6) = T^Tj^ is the cosine of the angle opened between the two vectors ki and /c2- The 

|fel||K2| 

need for a regularization ultraviolet cutoff to render the variance ()35|1 finite was already 
encountered in our discussion in Section III, case I. In next section, we will discuss this reg- 
ularization procedure in the context of a renormalization program of the physical parameter 
that measures the size of the mass anisotropics. 
Comparing with the general expression (jH) 

r2kuv d^C 

{AEyf = plj^ j^,m\)\Fv{0\\ (36) 
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where the geometric factor -FV(C) is given in expression Q, we can easily obtain the power 
spectrum of vacuum fluctuations: 

^(ICI) = f icos{9) - If, (37) 

in terms of the new variables = | (/I + C} and k2 = ^ (^jl — (j. The domain of integration 
in momentum space is defined by the condition d{Q = {fl E : km < ± C| < kuy}. 
This domain is the intersection region of two similar annulus with inner radius 2kiR and 
outer radius 2kuv and centered, respectiveley, in ±(f. 

It is easy to check that at C = 0, that is ki = k2, the power spectrum vanishes 

r{( = 0) = Um^_o P(C) = 0. (38) 

This result is a direct consequence of the fact that the quantum field is in its ground state, 
which is an eigenstate of the hamiltonian (j^Hj) and, therefore, the total energy E{Q) in the 
whole system Q does not fluctuate quantum mechanically, [AE^f = 0. This condition 
simply constrains the fluctuations of the density field to preserve the total energy (mass) in 

anisotropies (jH^ 



the system and, according to (j^^), it is sufficient to force the energy (mass 
to be scale invariant over cosmologically large covariant sub- volumes {AEvY oc S 

In fact, the scale invariance of anisotropies ()35|1 over cosmologically large comoving vol- 
umes was proved numerically in Additional numerical examples that beatufilly prove 
the claim were presented in In a first example, the variance of energy fluctuations of 
a free scalar field in 2 + 1 Minkowski space-time is calculated in the spatial region within 
a wrinkled surface. The interesting aspect of this example comes from the fact that the 
total volume contained within the surface does not change when the boundary surface wrin- 
kles, although the area of the boundary grows monotonically. The results presented there 
show that the variance of energy fluctuations does grow linearly with the growing boundary 
surface. In a second example, the variance of energy fluctuations is calculated within an 
annulus of inner radius ri and outer radius r2. When the radius ri gets larger, but still 
smaller than r2, the total volume of the annulus decreases. But the area of the boundary 
surface grows, and so also does the variance of energy fluctuations. 

When this mechanism of local rearrangement is discussed in the literature the attention 
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is focused on the quadratic coupling in the hamiltonian density ()33p of Fourier modes of 
the fields (pix) or 7c*{x) with close covariant momenta ki ~ /c2, which produces fluctuations 
of the density field T-C{x) with covariant wavelength A = 27r/|C| = 27T/\ki — that can be 
much longer than the horizon, even if both covariant wavelengths 27i/\ki\ and 27c/\k2\ are 
much shorter than the horizon. This is a very simple property of the coupling of Fourier 
modes that led to Y.B. Zeldovich, and many others after him, to notice many years ago 
that local fluctuations of the scalar field (jH^ with comoving wavelength shorter than the 
horizon can produce a non zero power spectrum ()37p 



^(ICI) - -oTTT^kuvlCl' (39) 



1 

Pl2(27r)2- 

over cosmologically short comoving momenta ( ^ H{teq). 

The interest in the literature on this aspect of the mechanism is motivated, nevertheless, 
by the erroneous assumption (fT^ . which wrongly asserts that mass (energy) anisotropies 
(I36|) over cosmologically large comoving volumes of size L ^ H^^{teq) are generated by the 
contribution dC, 'P(C)l-^v(C)P at ^ ~ L^^ <^ H{teq). This erroneous assumption (fT^ 
has misled the analysis of the mechanism of local rearrangement of matter as the origin of 
cosmological anisotropies since the early years when it was first studied. As we explained 
in Section III, case I), the contribution at ^ ~ to ()3(j|l is negligible compared to the 
ultraviolet contribution, d( '^(C)l-^v(C)P at C ~ kuv, in statistical systems globally 
constrained by the condition (AM^)^ = 0. This can be checked directly on ()35|) by noticing 
that the largest contribution comes from |A;i,2| ~ kuv when they are oppositely oriented, 
ki ~ — ^2, so that cos6 = —1, which contribute to the power spectrum V{( ^ ^/L) in the 
ultraviolet range, and not by the aligned modes ki ~ k2, which form an angle cos6 ~ +1 
and contribute to V{( ~ 1/L <^ H{teq))- When the dominant ultraviolet contribution is 
taken into account (fT3j) the resulting anisotropies over cosmologically large volumes are 
proven numerically to be scale invariant. 

The irrelevance of the quadratic dependence of the density field (jH^ in the fundamental 
fields (I32j) and (|34j) for generating scale invariant anisotropies over cosmologically large vol- 
umes out of fluctuations in Fourier modes of the latter with very short comoving wavelength 
is manifest in the next example. 

Let us consider the same theoretical framework laid in (jHTj) . (jH^ . but assume now that 
the matter density field is proportional to (fH^ the conjugate momentum field. Obviously, 
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the density field is still homogeneous and isotropic on average: 

(0|7r(f)|0) = 0. (40) 

But, as we know, this constrain cannot prevent the appearance of random quantum 
anisotropics in the spatial energy (mass) density distribution, because the ground state 
|0) is not an eigenstate of the restricted operators 

My (X I d^x 7r{x) + h.c. (41) 

if the volume of intergration V is not the whole 3D space Q. 

We find the variance of energy (mass) anisotropics in finite comoving volumes of cosmo- 
logical size to be 

{AMyf oc / f d'x dV(0|7r*(f)7r(y)|0) = ^ k |Fv(^)p, (42) 
Jv Jv JkiR [2tt)-^ 

given that {0\n*(x)n{y)\0) = J -^^k e*'^'^^"^). Now, the density field depends linearly in the 

fundamental fields, which are allowed only to randomly fluctuate in modes with comoving 

wavelength shorter than the horizon. We read from the expression above the power spectrum 

of the fluctuations in this example, V{k) ~ k for kj^ ^ k < kjjy and is zero otherwise. 

The ultraviolet cutoff kuv is necessary to render expression (02)) finite. Once we introduce 

the cutoff we immediately confirm, following previous discusssion in Section III, case II): 

(AMy)2 ~ (43) 

quantum anisotropics of the density field ()41|1 are scale invariant. This result was foreseeable 
on the ground (AMq)"^ = 0, simply expressing that also in this second example total energy 
(mass) is conserved and not allowed to fluctuate quantum mechanically. The need for an 
ultraviolet cutoff in ()42|) proves that these anisotropics are also mainly generated by the 
contribution from the modes with the shortest comoving wavelength. 

The two examples discussed in this section basically consist of random fluctuations that 
can displace the carriers of mass (energy), i.e. fluctuations of the fundamental fields ()32|l 
and ()34j) . only over causally connected distances. They both generate scale invariant mass 
(energy) anisotropics over cosmologically large comoving volumes. On the other hand, the 
two examples are defined by completely different density power spectra. The power spectrum 
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associated to ()35|) density fluctuations of the quadratic operator Tiv is 'P(fc) ~ fc^ over the 
whole range of cosmologically short momenta k < kuv-, with kuy > H{tf.q), while the power 
spectrum associated to (jl^ fluctuations of the momentum density operator My in the 
second example is V{k) ~ A; in the range of momenta within the horizon H{teq) ^ k < kuv- 
We already found this degeneracy in the power spectrum of scale invariant anisotropics in 
our previous analysis in Section III, case I). 

We close this section with a short comment on the cutoff procedures we have introduced 
above in this section. First, we realize that for the issues we have been discussing here the 
infrared cutoff km is irrelevant in both examples, as the contribution of the long wavelength 
modes to the total variance Q is sub-dominant. The dominant contribution comes in both 
examples from the short wavelength modes of the fundamental flelds. In the flrst exam- 
ple, fl35p . the dominant contribution is obtained from the quadratic coupling of oppositely 
oriented ultraviolet modes |A;i_2| ~ kuv, cos9 ~ —1. In the second example, the largest 
contribution to P2|l comes also from the ultraviolet modes of the momentum fleld. Both 
contributions are divergent and needs to be regularized with ultraviolet cutoffs. Now it is 
important to notice from ()19j) that the ultraviolet regularization procedure only affects the 
absolute size of the fluctuations, (AMy)^ ~ pi pS, but not the linear dependence on the 
area S of the surface that bounds the sub-volume. It is natural then to understand this 
procedure as an intermediate regularization step in the renormalization programme of the 
physical parameter p = y-Y . This comment introduces us to the issues addressed in next 
Section. 

VI. SCALE INVARIANT DENSITY ANISOTROPIES IN QUANTUM FIELD 

THEORIES. 

The scalar quantum fleld theory (j^T|) of last section is obvioulsy only a toy-model de- 
signed to confront some widespread, but wrong, beliefs about the dynamics of primordial 
structure formation. The model has also been a satisfactory bench to discuss some impor- 
tant concepts that shall be very relevant in the development of a complete quantum fleld 
theory of primordial cosmological density perturbations. 

Within the setup of a QFT the energy (mass) density fleld at the time of equality ((H) is 
described by a density operator 
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^ h ^'""'^ (44) 

analogous to the density operator in p3|l . acting on the quantum state of the system |0). 
This state is required to be invariant under spatial translations and rotations. Homogeneity 
and isotropy demand (0|/i(x)|0) = (0|/;o|0), which implies 

(0|%|0) = 0, (45) 

and, therefore, the average value of M(y) = Jyd^^ h{x), the energy (mass) contained 
within a spatial macroscopic sub-volume V, is proportional to the volume of the region, 
(0|M(l^)|0) = V X (0|/io|0). The magnitude M{V) is said to be an extensive properties of 
the system. We assume here, as we did in Section III, that {ho) = (0|/io|0) 7^ 0. 

Condition ()45|) does not prevent, nevertheless, the appearance of random spatial energy 
(mass) density anisotropics due to quantum fluctuations, with variance 

{AEy)' ^ (0| [M{V) - {M{V))]' |0) = {h^fj ^ VmFvm', (46) 
where the power spectrum V{k) is defined, as usual, by 

{0\5l 5^,J0) = {27rfV{h)6'{h - h). (47) 
In a wider context, (jl^ can be the density field operator of any global Noether charge 

n = [ h{x) = f k). (48) 
Jn Jn 

The physical observable M{V) describes then the net value of that charge within a finite 
macroscopic comoving sub-volume of cosmological size. Condition implies that charge 
H is, on average, homogeneously and isotropically distributed in space, while reminds 
us that quantum fluctuations can, nevertheless, produce random anisotropics in its spatial 
distribution. 

According to the central limit theorem, for a complete statistical description of random 
anisotropics in the spatial distribution of the extensive charge is enough to provide the value 
of their variance ()46|) over comoving volumes of arbitrary cosmological size. The calculation 
of the variance of random spatial anisotropics in the distribution of a Noether charge is a 
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well-posed physical question that should get an unambiguous, finite answer in any properly 
formulated theoretical framework. Of course, the theoretical calculation must produce the 
clear result [AEvY = over any macroscopically large comoving volume if and only if the 
system does not generate random anisotropics. 

Hence, it is intriguing to notice that within the setup that we are considering we are 
allowed to modify the definition of the charge density field operator ()44j). for example by 
adding a total three-divergence, without changing the definition of the global charge ()48p . 
Such 'symmetry' transformations of the charge density field ()44|) must preserve the definition 
of the zero-mode ho and they must also respect the expectation values ()45j) of the non-zero 
Fourier modes 7^ in order to preserve the average homogeneity and isotropy of the 
system, but they can alter the power spectrum V{k) at A; 7^ 0, defined in (jUj). Therefore, 
they can alter the variance of quantum mechanically generated spatial charge anisotropics in 
finite macroscopic sub- volumes, defined in ()46|). This degeneracy implies, in principle, that 
these two concepts (j47jl and ()46|1 cannot be uniquely defined in the considered theoretical 
framework and raises questions on the capability of the setup to properly describe quantum 
mechanically generated spatial anisotropics. In the paragraphs below we investigate and 
clarify this point. 

We need first to identify the invariant features of the power spectrum (jlTj) and variance 
PUj) under the 'symmetry' transformations of the density field (jl^ : as the zero mode ho is 
not modified by the permitted transformations of the density field operator, its statistical 
momenta (0| {ho — {ho)y |0), p = 1,2,3, are invariant features. The first moment p = 1 
is zero, by definition. The second and higher momenta, on the other hand, are zero only if 
the quantum state |0) is an eigenstate of the global charge (^Hj). In particular, this implies 
that energy (mass) is globally conserved, 

(A^n)' = 0, (49) 

by quantum random fiuctuations of the density field (03)). 

This condition ()49|1 is naturally satisfied if |0) is the ground state of the system and 
is a conserved Noether charge. According to our result ()24|1 in section IV such constrain 
implies that quantum mechanically generated anisotropics (I4()j) in the spatial distribution of 
charge 7i must be scale invariant Q over comoving volumes of cosmological size. Thus, the 
linear dependence of the variance of primordial cosmological energy (mass) anisotropics on 
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the area of the surface that bounds the considered comoving volume, 



Ln {AEv.f = Ln {AEy^f + Ln 



( 



5i 

^2 



) 



(50) 



where Vi and V2 are arbitrary macroscopic sub- volumes and 5*1, S2 are, respectively, the 
areas of the surfaces that bounds each of them, is an invariant feature under 'symmetry' 
transformations of the density field operator in the ground state of the system [ij]. This 
is nothing but an example of the basic feature that we already noticed in our analysis of 
Section III, case I), that power spectra with different spectral index n > 1 cannot be distin- 
guished through the energy (mass) density anisotropics they produce over macroscopically 
large comoving volumes: 'symmetry' tranformations of the density field ()44j) simply modify 
the power spectrum ()47|1 of spatial anisotropics 'P{k) ^{k) + AV{k) by a regular term 
AV{k) = o{k), which does not alter the macroscopic scale invariance Q of the anisotropics. 

The transformations of the density field can modify, on the other hand the size of the 
spatial charge anisotropics over macroscopically large spatial sub- volumes, parameterized by 
the factor of proportionality fi, 



This parameter /i is not an invariant feature under 'symmetry' transformations in the def- 
inition of the charge density and we should conclude that it describes a physical property 
that is not properly defined in this setup. We explain below the origin of this problem and 
argue, using usual statistical/QFT renormalization concepts, that this parameter fi must be 
considered as an additional free parameter of the theory, like masses or coupling constants, 
which fixes the size of quantum mechanical random spatial charge anisotropics in the ground 
state of the system. 

In Section III, case I), we learned that scale invariant energy (mass) anisotropics over 
macroscopically large volumes are dominated by the contribution of the ultraviolet modes 
to the integral expression ()17j) . To render this expression finite we introduced by hand an 
ultraviolet cutoff, which fixes the size of the mass (energy) anisotropics in cosmologically 
large comoving volumes. We can now understand this step in Section III as a regularization 
procedure of the infinite parameter /i. The 'symmetry' transformations of the density field 
can be reabsorbed within this infinite parameter. Or, in a different point of view. 




(51) 



32 



the infinite parameter n can be regularized by adding an appropriate counter-term in the 
density field. This procedure do not modify the finite relationship (j^Uj) . which is a physical 
prediction of the setup. In conclusion, scale invariance of vacuum anisotropics, expressed 
through the finite relationship must be respected through the renormalization program. 
The price we pay is the normal price in QFT/statistical mechanics: we cannot predict from 
first principles within this setup the value of the 'dressed' parameter /i, in the same way 
that we cannot predict the physical value of in QED. 

Only after fixing, normally through measurement, the variance {AEy^y = A^, of ran- 
dom anisotropics over a macroscopic volume of reference Vi , we have a meaningful physical 
prediction 



for the variance of mass (energy) anisotropics over any other comoving volume V2 of cos- 
mological size. This relationship is equivalent to the renormalization group equation that 
predicts physical values for the running of the coupling constant aem{p'^) at arbitrary energy 
scales only after an initial value aem(Po) a certain renormalization energy scale has been 
fixed by experiment. 

Please notice that the incapability of the QFT setup to predict the 'dressed' value of 
the physical parameter /i due to common ultraviolet divergences and the impossibility to 
single out in the setup the local definition of the energy (mass) density operator, reaches 
also the widely accepted predictions of the infiationary scenario. We can freely add a three- 
divergence to the hamiltonian density, which will dress the scale invariant anisotropics pre- 
dicted by infiation. As we have said, the additional terms respect the scale invariance of 
the anisotropics but make the parameter fi in (jHTjl formally infinite. Let us say it in 
different words: under mistaken assumption (fTT|) . upon which the calculations of the cosmo- 
logical anisotropics in inflationary cosmology are carried out, there are no divergences in the 
theoretical calculation of the anisotropics, but just because the inflnites appear only after 
noticing the degeneracy (|15|) in the power spectrum of scale invariant mass (energy) cosmo- 
logical anisotropics. Ordinary predictions of inflationary cosmology of primordial structure 
does not suffer from the problem of inflnite expressions, only because in these calculations 
the inflnities are thrown away, using approximation (jlll) . without justiflcation. A correct, 
not based on mistaken assumption complete calculation of inflationary predictions will 




(52) 
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also lead to finite relationship ()5U|). but only after a regularization procedure of the infinite 
expression for the parameter fi. 

It is very interesting to notice that the infinite contribution from ultraviolet modes of the 
density field to the physical parameter /i in ()5H1 is, according to (j23), directly associated 
to the fact that total mass (energy) is globally conserved (AEq)"^ = 0. Scale invariant 
anisotropics of extensive magnitudes defined in the bulk of an spatial volume V are associated 
to the surface that bounds the sub-volume. If these anisotropies can be holographically 
described by a theory defined on the surface, as suggested in interestingly find 

a connection between the infinite, renormalizable parameter /i of the holographic theory 
defined on the closed surface, and the globally conserved extensive magnitude defined in the 
interior and exterior volumes to that surface. In quantum Hall effect of strongly correlated 
electrons confined in 2D spatial dimensions by a perpendicular magnetic field in condensed 
matter physics it is known that is possible to obtain an holographic ID description of the 
edge excitations over the boundary border 

The problem of determining the size n of scale invariant energy (mass) anisotropies 1)5111 is 
not very different from the problem of determining the cosmological constant in the context 
of relativistic QFTs. The average value {M{V)) = V x (/iq) is, in general, ultraviolet 
divergent. We cannot predict within the QFT setup the value of A = (/iq), but once the 
theory is regularized we find that (-Evi) = ^(-^^2)? the energy (mass) contained within a 
macroscopic sub-volume is an extensive magnitude. This is a finite relationship that must 
be respected whatever the unknown mechanism that renormalizes the value of A could 
be. In the same way we cannot predict the physical value of the parameter /i, but we 
can predict the scale invariance ()50|) of the variance of spatial charge anisotropies. Both of 
these problems become physically meaningful only when the energy (mass) density field gets 
coupled to gravity and, therefore, they can only get a final answer in the more fundamental 
setup that includes the quantum description of both matter components and dynamics of 
the space-time metric. 
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VII. SCALE INVARIANT ANISOTROPIES IN THE GAUGE INVARIANT 
FORMALISM OF LINEARIZED DENSITY PERTURBATIONS. 

Spatial anisotropies in the energy-momentum density field necessarily perturbate the ho- 
mogeneous and isotropic FRW expanding space-time background, which in turn can feed 
back the dynamics of the density anisotropies. A complete analysis of the locally Lorentz 
invariant equations of general relativity that describe the coupled dynamics of density and 
metric perturbations is technically very demanding even at the classical level because it must 
carefully take into account unphysical gauge degrees of freedom. Fortunately, a first aprox- 
imation of the coupled dynamics at linear order in the perturbations has been formulated 
in a very simple ^auge invariant framework. A compiling detailed report of this approach 
can be found in The linear approach is justified in the literature on the argument that 
cosmological primordial mass (energy) perturbations are tiny, 5p/ p < 10~^. In the previous 
section we have found, nevertheless, that in the standard statistical setup the physical pa- 
rameter /i, which measures the size of the fluctuations, 5p oc y/p S, is infinite and needs to 
be renormalized. This brings us to face another usual procedure in perturbative QFT: the 
perturbative expansion must be carried out in powers of the dressed physical value of the 
parameter 6p/ p, which is really tiny, and not in terms of the infinite, or very large, undressed 
parameter. We will not further discuss here how to carry a renormalization programme. In 
this section we only explore within the gauge invariant framework at linear order the con- 
cepts introduced in the previous sections. This section VII, like Section VI, report only first 
results of research in progress. 

First, let us very briefly review the basics of the linearized gauge invariant formalism, 
which start from the usual Einstein-Hilbert action for the coupled dynamics of the space-time 
metric and energy-momentum matter tensors. Both tensor field are then decomposed as a 
sum of two terms: a first homogeneous and isotropic FRW background term plus a second 
perturbative term describing the anisotropies. Using the equations of motion the action 
is rewritten in linear aproximation keeping up to second order in the perturbative terms. 
After some straightforward, but lengthty, calculations the general action of coupled density 
and metric perturbations during the radiation dominated epoch of a FRW background is 
written in comoving coordinates in terms of a single gauge invariant scalar field potential 
v{x), whose dynamics is described by the hackneyed hamiltonian of a free scalar field in fiat 



35 



Minkowski space-time 



n = f n{x)= f d^x(n{x)n{x) + lvv{x)Vv{x) 




(53) 



Field operator v{x) and its conjugate momentum Tr{y) obey canonical commutation relations 
[v{x),-iT{y)] = i 5^{x-y). 

During the radiation dominated stage of the standard cosmology there is not cosmological 
particle production and, therefore, the gauge invariant scalar field remains in its ground state 
|0) when the temperature of the universe cools down to the onset of matter radiation equality. 

Furthermore, every physical observable can be expressed in this formalism in terms of 
this pair of conjugate fields, thus eliminating the problem of spurious gauge degrees of 
freedom. In particular, the Fourier modes of the energy (mass) density field operator 
with comoving wavelength shorter that the horizon are proportional in this formalism to the 
Fourier modes of the momentum density operator: 



where Teq is the cosmic temperature at the time of equality and diJ^x) is an undefined 
arbitrary three divergence according to our discussion in the previous section. Formally, 
this density field is the second of the two examples we discussed in the toy-model of Section 
V. We found there that in the ground state of the system |0) its spatial anisotropics are 
scale invariant (021) • 

We consider important to remark that the density field operator in the linearized gauge 
invariant formalism is proportional to the momentum density (jH^ only over causally con- 
nected distances. Over super-horizon distances they do not necessarily coincide and the right 
hand side of ()54|) gets correction terms. This is not really important to us, as we know that 
the largest contribution to scale invariant cosmological anisotropics, globally constrained to 
preserve the total energy (mass) of the system, comes from fluctuations in under-horizon 
modes of the density field (fTTH) . These modes are described by (jK^ . 

The unsolved situation regarding the theoretical prediction of the parameter in (15111 
has not improved in this formalism of linearized density anisotropics coupled to metric 
perturbations. A direct estimation of the variance of energy (mass) anisotropics from ()54|) 
produces 




1 



1 



(54) 
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{^Eyf ^ {h^Y^S. (55) 

eg 



Over a sphere of comoving radius L we then obtain 



5pV (AE- 



V, 



\2 



\pjy {hyV^ 477 (Te,L)4- ^^^^ 

A rapid estimation for L ~ 10^H~g^ ~ 10^^ gives for the un(iressec? dimensionless parameter 
{^^y ~ (^) ~ 10^°°. Obviously, this is much larger than the observed physical 

dressed value (^)^^ ~ 10~^. The renormalization can be carried out through counter-terms 
introduced by the three-divergence in the definition of the density operator ()54|) . 

For a deeper treatment of this problem is necessary to go beyond the linear order of the 
density and metric perturbations. We guess that the final quantum theory of gravity has 
a say in fixing the size of scale invariant primordial anisotropics. Dynamics of structure 
formation could thus offer a bench for phenomenologically testing that final theory at the 
low energies of matter-radiation equality. 



VIII. RENORMALIZATION EQUATIONS OF PHYSICAL DENSITY 

ANISOTROPIES. 

We wish to explore in this section how the renormalized (dressed) mass (energy) 
anisotropics depend on the resolution scale at which the physical density field is probed 
and obtain, within the theoretical setup laid in section VI, the renormalization group equa- 
tions that describe this dependence. 

Let us say that Pphy{x; Iq) is the physical density field probed with arbitrary length reso- 
lution Iq = Mj"^. When the system is probed at blunter length resolution li = Mf ^ > Iq (or, 
in energy scale Mi < Mq), the observable field 

Pphy{x] k) = / d^y pphy{m h) W{y- x; k, Iq) (57) 

looks further smoothed because its Fourier modes with wavelength A < /i shorter than the 
size of the new probe are cut off by an appropriate convolution kernel (window function) 
W{y — x; li, Iq) extended around the point of observation x over a certain domain of typical 
size li and normalized to 
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/ d'yW{y-x;h,l^) = l. (58) 

Expression (|57j) simply means that measuring a local observable at any point x G with a 
probe of length li does in fact return an average weighted value of the field over a vicinity of 
the size of the probe around the tested point. This is typically the shortest scale that that 
probe can test. 

Obviously, the physical density field (jFTj) resolved to lengths of the order li cannot depend 
on the arbitrary shorter scale Iq used in the definition above. This scale Iq plays here a 
role similar to the arbitrary renormalization scale in renormalization group equations. The 
condition of independence of the physical observable field ()57j) resolved at length li on the 
renormalization scale k) fixes the functional dependence of the convolution kernel 

/ d''zW{y-z;l,lo)W{z-x;h,l) = W{y-x;h,k>). (59) 
Jn 

The simplest solution to this convolution equation is a Dirac delta function W{x — y] la, 4) = 
6^{x — y) for whatever two resolution length scales 4, 4 we choose to compare, but it must be 
understood that this choice is rather unphysical as it assumes that we can solve the density 
field with infinitely sharp length resolution. A more realistic solution to is a normalized 
gaussian kernel: 

W{x-y-,L,h)=ArExp[-L^]=ArExp[- ^,'_^^ \x-y\% (60) 

where 

Naturally, when Ma -C Mb the convolution kernel is W{x — y; la, k) — N'Exp[—M'^\x — ^p]. 
On the other hand, when we compare the observed physical density fields at two very close 
resolutions scales 4 ~ 4, the convolution kernel is roughly W{x — y; la, 4) — S^{x — y), as it 
should be. 

It is also useful to compare in momentum space the physical density fields resolved at 
different length scales. Expanding Pphy{x; li) and Pphy{x; Iq) in Fourier modes, as in eq. ([T]), 

5g[Mi] e-''-', (62) 
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6^[Mo] e-'^-"^ (63) 



and then introducing the second expansion (jUHj) in eq. ()57|) 



Pphy{x; h) = 1^ \po + PoJ ^ S^,[Mo] e-^'-y ] W{y- x; h, 4,) 



W) 

Thus, we obtain by comparison with (j62j) the relationship 



6^[M^] = W^[M^, Mo]%[Mo] = e 6^[Mo], (64) 



where 



W^[Mi, Mo] = ^ d'y e-'^<y-^^ W{y- x; h, /q) = e ' k '^i ' , (65) 

is the Fourier transform of the convolution kernel. We find ()64p that, as a result of changing 
the scale at which we resolve the density field, the stochastic modes get renormalized 
by the window function ()65|1 . Therefore, the power spectrum associated to these stochastic 
modes, defined by (27r)'^'P(/c) b'^{k\ — k^) = (51 5^^), gets renormalized as 



V,hyikuM,) = \W^,^[M,,Mo]\'V,hyik,;Mo) = e ' ^'o ^'t V.^yih; Mo). (66) 
If we keep the arbitrary scale Iq fixed for comparison and let the resolution scale 4 = M{'^ ^ 



/o to run, Vphyiki, Mi) ^ e ^ "i Vphy{ki; Mq), we clearly see that the resolution scale simply 
introduces a running physical ultraviolet cutoff Mi, beyond which the power spectrum of 
the density fluctuations is exponentially suppressed. 

With these tools at hand we can obtain the variance (0)) of physical density anisotropies 
(AMy[M])2 = J ^Vphy{k; M)\Fv{k)\'^ as a function of the scale / = M'^ at which the 
density field (jFTj) is resolved. The important feature to be noticed is that, for whatever two 
resolution scales Mi,Mo that we choose to compare, the normalization condition (jKH|) fixes 
Wk=o[Mi, Mq] = 1 at the origin = in momentum space, see eq. (jHHjl . Hence, the power 
spectrum at the origin is invariant Vphyik = 0;Mi) = Wk=o[Mi, Mq] Vphy{k = 0; Mq) = 
'Pphyik = 0; Mo) under running of the resolution scale M. In particular. 
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Pp,,(A; = 0;Mi) = ^ ^^.^(A; = 0; Mq) = 0. (67) 

According to eq. (j22|l . = 0) = is a necessary and sufficient condition for the variance of 
mass (energy) anisotropies in a finite sub- volume V to be scale invariant. Let us remind that 
dV{k)/dk\k=o = is necessarily fulffiled ifV{k = 0) = 0, up to some technical considerations. 
Therefore, we must conclude that scale invariance of mass (energy) anisotropies in finite 
macroscopic sub- volumes is an invariant feature of the statistical density anisotropies under 
transformations of the scale / at which we resolve the density field (P), as long as this scale 
is typically shorter than the size L of the considered sub-volume / ^ L. 

In eq. fl24j) we related the scale invariance of mass (energy) anisotropies in macroscopic 
sub-volumes to statistical fluctuations of the density field ()57|1 constrained to conserve the 
total mass (energy), (AMq)^ = 0. Equation (jU7j) above simply expresses that this global 
constrain should remain unchanged when we change the resolution scale at which we test 
the density field. 

It is interesting at this point to bring back our analysis of Section III, case I. There we 
found a general approximate analytic expression for the variance of scale invariant mass 
(energy) anisotropies in finite macroscopic sub-volumes of size L, after introducing by hand 
a regularization scale kc- 

{AMvr - Apl fj' dikL)^ = Apl dk^, (68) 

which we wrote in ^ as (AM^)^ = pi pS, with p = A Jq'' dk^^. We found that the 
parameter p is formally infinite when we take the regularization scale kc to be infinitely 
large and noticed in Section VI that this divergent parameter can be renormalized to its 
physical value following an appropriate procedure. It is also worth to comment that in our 
discussion in Section III the regularization scale kc was introduced as a sharp ultraviolet 
cutoff in the integral above, but we remarked there and in Fig. 1 that similar results are 
obtained when we use, instead, exponential or polynomial cutoff procedures. 

We now find exactly the same formal expression for the renormalized (physical) variance 
of mass (energy) anisotropies in finite sub- volumes V of typical comoving size L when the 
fluctuating density field is probed to a length scale / = M~^. 
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(AMyYlM] ~ Apl L'f'dikL)'^^^ = Apl [ dk^^^. (69) 

In this expression the upper hmit in the integral is not a sharp ultraviolet cutoff, but it 
simply means that the resolution scale introduces a physical exponential cutoff at 
k > M. This expression is obvioulsy finite. It describes a physically observable magnitude 
(a dressed magnitude in the common vocabulary of renormalization group equations). It is 
also clear from this expression that the effect of changing the resolution at which we probe 
the mass (energy) anisotropics can be reabsorbed in the renormalizable running parameter 
fiphy{M) = 4 dk^^^^^^, which measures the typical size of the scale invariant anisotropics 
at that resolution scale [22], (AMy)2[M] ~ pi pphy{M)S{V), such that 

= -M' ^Im = -4 Vp,y{M). (70) 

This equation is a renormalization equation in the most common sense, but it deserves 
an explanation. We have discussed in Section III that all power spectra with index n > 1 
are physically indistinguishable through the scale invariant anisotropics (AMy)^[M] they 
produce over finite, but macroscopically large, volumes V because, we said, the different 
indices n can be reabsorbed in the renormalized parameter Pphy The analysis in this section 
corroborates this conclussion and it offers us a deeper insight. At a given resolution scale 
all the information contained in any power spectra with index n > 1 is hidden in the 
renormalized parameter Pphy, which measures the overall size of the scale invariant physical 
anisotropics. As we cannot predict the physical value of this parameter in the theoretical 
setup that we have laid, different power spectra with n > 1 are physically indistingishable 
at this point. On the other hand, the spectral index n can be physically observable through 
the renormalization equation (f7n|) . that is, through the change in the overall size of the 
scale invariant mass (energy) anisotropics when we let the resolution scale M to run. We 
must conclude that, even though all power spectra with index n > 1 produce scale invariant 
mass (energy) anisotropics over macroscopically large volumes, different spectral indices are 
physically distinguishable if we can reliably measure the change dpphy/dl in the overall size 
of the anisotropics over a given macroscopic volume when we change the resolution scale 
/. This is the reason why in this section we have labelled the power spectrum Vphyik) as a 
physically observable magnitude. 
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An additional note. In the preceding paragraphs we have discussed how the resolution 
scale at which the random density field (jFFj) is probed affects the theoretical predictions on 
physically observable mass (energy) anisotropics. The analysis led us to the renormalization 
equation ()7()|) . Quite naturally we found that the resolution scale introduces an exponential 
physical cutoff that suppresses all fluctuating modes with comoving wavelength shorter than 
the chosen scale of resolution. If this comoving scale is conveniently shorter than the co- 
moving size of cosmological volumes the observed mass (energy) anisotropics are necessarily 
scale invariant. 

Before closing the discussion on this issue we want to comment on a different approach 
about how to manage the dependence of physical observables on the flnite resolution scale 
of any real probe. This second approach is widely discussed in the literature on primordial 
anisotropics, but it is somehow artiflcial. 

Let Pphy{x) be the physical density fleld resolved at a certain scale Iq. The mass in a 
macroscopic sub- volume V of typical size Iq (which we assume for simplicity centered at 
the origin x = 0) is deflned in this second approach as MgaussiV) = In d'^^Pphy{x)e~^^^'^^^'^ , 
where A/" normalizes the gaussian weight, instead of the deflnition M{V) = Jy d^xpphy{x) we 
introduced in Section II and have used throghout the paper. The new operator MgaussiV) is 
commonly referred in the literature as the mass function with a normalized gaussian window, 
while M{V) is commonly referred as a top-hat mass function. 

In Section VIII we have discussed how to properly accomodate a flnite resolution scale 
within the standard definition of a top-hat mass function. The use of a mass function with a 
gaussian window Mgauss{V) follows a different, and more artificial, approach: it averages the 
density field over the whole 3D space using a weighting gaussian function Af e~^^/^^ with 
a typical opening width of the size L of the considered volume V. According to eq. (j37j) 
the gaussian mass function MgaussiV) oc Vpphyix = 0;L) is nothing but the average mass 
in volumes of size L when the physical density field is probed with a resolution scale of the 
same size of the volume V of integration, / ~ L. The aim behind the definition Mgauss{V) is 
to introduce in the theoretical expression for the variance of mass anisotropics in a sub- 
volume of typical size L an ultraviolet cutoff in momentum space of the order ~ 1/L and 
artificially justify the estimation It is important to notice that within this gaussian 

approach the volume V is not treated as a macroscopic volume, but as a volume of the 
minimal size that the probe can test. And we must be aware that a theory endowed with 
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an ultraviolet cutoff in real space of size / can lead to contrieved predictions when applied 
to physical sub-systems whose typical size is L < /, in the same sense that the effective 
Fermi four-point vertex fails to reproduce the observational features of weak interactions at 
energies close or higher than the mass of the W and Z bosons. In both cases, we might be 
trying to apply the theory beyond the scales it can reach. In the case of a power spectrum 
with index n = 1 the result of using the gaussian mass function Mgauss{V) does not modify 
the scale invariance ()5U|) . but it artificially render the divergent parameter /i in ()51|) finite. 
In the case of a power spectrum with index n > 1 the use of a gaussian mass function 
artificially produce the wrong results (fT^ . instead of the scale invariant (|Hn|l . 

IX. DISCUSSION. 

We have revised in this paper the standard theoretical framework commonly used to 
discuss the origin of scale invariant cosmological mass (energy) anisotropies in the homo- 
geneous and isotropic FRW background at the time of matter radiation equality to point 
out a crucial fault (fT^ in the argument that led to wrongly assert the so-called origin of 
structures problem of standard cosmology. We have shown that the correct evaluation (fT^ 
offers a natural and elegant explanation of the scale invariance of the primordial cosmolog- 
ical anisotropies as a consequence of a conservation law ()24j) of the total mass (energy) of 
the fluctuating density fleld. Surprisingly, the correct estimation (fTH|l that we present here 
is well-known since long ago in statistical mechanics and condensed matter physics, but was 
never properly noticed before in cosmology. 

According to our analysis primordial cosmological mass (energy) anisotropies in the, 
otherwise, homogeneous, isotropic and flat FRW universe at the time of matter radiation 
equality could happen to be simply quantum vacuum fluctuations of the density fleld gen- 
erated at the same instant of equality with comoving wavelength shorter than the causal 
horizon at that time. In other words, the mass (energy) primordial anisotropies over cosmo- 
logically large volumes would happen randomly as a result of local rearrangement of matter 
at the instant of equality through the surface that bounds the considered spatial regions, 
thus explaining why their variance is proportional to the area of the boundary surface. 

Obviously, this explanation of the origin of primordial cosmological anisotropies at the 
time of matter radiation equality in the context of standard cosmology does not require 
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any previous epoch of inflationary expansion. In the context of this alternative scenario 
the fluctuating density field is naturally required to be in its ground state, thus avoiding 
also the problematic issues on the exceptional initial conditions of the universe that can 
produce a period of inflationary expansion Furthermore, the physics involved in the 

generation of the anisotropies shall probably be physics at the scale Teg ~ 1 eV of equality, 
instead of the very high energy scales T > 10^^ GeV summoned in inflationary mechanism. 
Yet, this alternative mechanism to the origin of large scale cosmo logical structures needs a 
further elaboration of the calculations presented in the last section VII on the interplay of 
the surface mass (energy) anisotropies with the theory of gravitational structure formation 
beyond the linear approximation. 

Finally, we must remark that the correct estimation which lies at the foundation 
of all the arguments in this paper, implies a continuous degeneracy in the power spectrum 
(H)) of statistical systems that produce scale invariant mass (energy) anisotropies, which we 
expressed in our discussion through the impossibility to distinguish spectra with indexes 
n > 1 over cosmologically large comoving volumes. The degeneracy is associated to the 
possibility of transforming the charge density field operator of a globally conserved Noether 
charge without actually modifying the charge and it has far reaching consequences. Such 
transformations do not modify the scale invariance of the spatial quantum anisotropies in the 
distribution of the charge, but they can actually change the formal expression of the infinite 
parameter that measures the size of these anisotropies. Or, from a different perspective these 
transformations can be reabsorved in the definition of an infinite renormalizable parameter 
of the theory. From this point of view we understand the scale invariance statement (jKn|l 
in the sense of a renormalization group equation and conclude that the physical dressed 
parameter, i.e. the absolute size of the anisotropies, cannot be predicted in the available 
setup of quantum field theories. In the same sense that we cannot predict the actual value 
of coupling constants or masses, or the cosmological constant. 

The picture of the origin of primordial cosmological fluctuations as vacuum fluctuations 
at the same instant of matter radiation equality would not be complete without a discussion 
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, advanced that in the context of 
20l fluctuations of local densities 



of how vacuum fluctuations can decohere. Halliwell, in 
the decoherent histories approach to quantum mechanics 
(momentum, energy) are more prone to decohere. More recently we have developed a new 
formalism j2l| to address the same issue: in this formalism the ground state of a bosonic or 
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fermionic system is described as a linear combination of randomly distributed pseudoclassical 
incoherent paths, which allows a description of the fluctuations of collective operators like 
the energy in a sub- volume in terms of classical stochastic concepts. 

X. ACKNOWLEDGMENTS. 

I am glad to thank J. Oaknin for his suggestions. I am thankful to A. Yarom and 
R. Brustein for discussions on a related project. I am also thankful to M. Joyce and D. Chung 
for their comments on an earlier version of this paper. This work was partly supported by 
the National Science and Engineering Research Council of Canada. I am grateful to Ana 
Oaknin, my teacher, tutor and friend, in memoriam. 



[1] E.R. Harrison, Phys. Rev. D 1, 2726 (1970); Ya.B. Zeldovich, Mon. Not. R. Acad. Soc. 160, 
1 (1972). 

[2] C. L. Bennett et al, Astrophys. J. 436, 423 (1994) | arXiv:astro-ph/9401012| . 

[3] D. N. Spergel et aL, |arXiv:astro-ph/Q302' 209 ; E. Komatsu et al, ^arXiv:astro-ph/Q302'223 

[4] C. B. Netterfield et al. [Boomerang Collaboration], Astrophys. J. 571, 604 (2002) 

[arXiv :astro-ph/ 0104460 ; C. Pryke, N. W. Halverson, E. M. Leitch, J. Kovac, J. E. Carlstrom, 

W. L. Holzapfel and M. Dragovan, Astrophys. J. 568, 46 (2002) a rXiv:astro-ph/0104490, ; 

R. Stompor et al, Astrophys. J. 561, L7 (2001) arXiv:astro-ph/0105062| . 
[5] M. S. Turner, "The new cosmology: Mid-term report card for inflation," 

larXiv: astro-ph/02 12"28T 
[6] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992); 
[7] P.J.E. Peebles, The large scale structure of the universe (Princeton University Press, NJ, 

1980); T. Padmanabhan, Structure formation in the universe (Cambridge University Press, 

Cambridge, England, 1993); S. D. White, Lectures given at Les Houches, August 1993, 

|arXiv:astro-ph/9410043| 
[8] B.I. Halperin, Phys. Rev. B25 (1982), 2185. 

[9] Ya. Zeldovich, Adv. Astron. Ap. 3, 241 (1965); Ya. Zeldovich and I. Novikov, Relativistic 
Astrophysics, Vol. 2, (Univ. Chicago Press, 1983). 



45 



[10] S. Hollands and R. M. Wald, Gen. Rel. Grav. 34, 2043 (2002); L. Kofman, A. Linde and 



V. Mukhanov, JHEP 0210, 057 (2002); S. Hollands and R. M. Wald, |aIXh^ihep-tli/0210001, 
[11] J. Beck, Acta Mathematica 159, 1-878282 (1987) 

[12] T. Bringmann, C. Kiefer and D. Polarski, Phys. Rev. D 65, 024008 (2002). 

[13] A. Gabrielli, M. Joyce and F. Sylos Labini, Phys. Rev. D 65, 083523 (2002) 
|arXiv:astro-ph/011045l| ; A. Gabrielli, B. Jancovici, M. Joyce, J. L. Lebowitz, L. Pietronero 
and F. Sylos Labini, Phys. Rev. D 67, 043506 (2003) arXiv:astro-ph/0210033] ; 

[14] D. H. Oaknin, |aIXiv:hep-th/0305068 l 

[15] B. Carr and J. Silk, Astrophysics Journal, 268, 1 (1983); J. Robinson and B. D. Wandelt, 



Phys. Rev. D 53, 618 (1996) |arXiv:astro-ph/9507043| . 
[16] R. Brustein, D. Eichler, S. Foffa and D. H. Oaknin, Phys. Rev. D 65, 105013 (2002) 

|arXiv:hep-th700 09063 . 
[17] R. Brustein, D.H. Oaknin and A. Yarom, |hep-ph/031009ij 
[18] R. Brustein and A. Yarom, |arXiv:hep-th/0302186| 

[19] J. J. Halliwell, Phys. Rev. D 58, 105015 (1998) |arXiv:quant-p h/9805062|. 

[20] M. Gell-Mann and J. B. Hartle, Phys. Rev. D 47, 3345 (1993) |arXiv:gr- qc/9210010|; 

R.B. Griffiths, Phys. Rev. Lett. 70, 2201 (1993); R. Omnes, Rev. Mod. Phys. 64, 339 (1992). 
[21] R. Brustein and D. H. Oaknin, Phys. Rev. D 67, 025010 (2003) arXiv:hep-th/0207251, ; 

D. H. Oaknin, Phys. Rev. D 67, 105016 (2003) |arXiv:hep-th/0304078| . 
[22] We are assuming at this point that the considered sub-volume V is bounded by a smooth 

enough surface, such that its area S{V) does not change with the resolution scale at which we 

probe the 3D space. 




FIG. 1: The points on the plot represent the numerically evaluated normalized ratio ^ ^/"^ of the 
variance of mass (energy) fluctuations in spherical 3D volumes of comoving radius L to the area of 
its surface S = AttL'^. Two different power-law spectra V{k) = ak^ with positive spectral indices 
are considered: n = 1 (dashed lines) and n = 4 (straight lines). The normalized data clearly show 
that in both cases the ratio remains constant over a very large range (100 < — > 10000) of values of 
the comoving radius L, measured in arbitrary units of length. The cutoff scale has been chosen 
arbitrarily at k^^ ~ 0.1 units of length, so that the dimenssionless parameter kcL ^ 1 ranges 
from 10^ to 10^. Two different setups to cutoff the power spectrum in the ultraviolet modes has 
been considered in this graph: an exponential cutoff P(A;)Exp(— fc/fcc) (squares); and a polynomial 
cutoff :P(fe) (circles). 



